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Abstract —Motivated by the fact that full diversity order is 
achieved using the ’’best-relay” selection technique, we consider 
opportunistic amplify-and-forward and decode-and-forward re¬ 
laying systems. We focus on the outage probability of such a 
systems and then derive closed-form expressions for the outage 
probability of these systems over independent but non-identical 
imperfect Log-normal fading channels. We consider the error 
of channel estimation as a Gaussian random variable. As a 
result the estimated channels distribution are not Log-normal 
either as would be in the case of the Rayleigh fading channels. 
This is exactly the reason why our simulation results do not 
exactly matched with analytical results. However, this difference 
is negligible for a wide variety of situations. 

I. Introduction 

Cooperative communication increases the performance 
quality of communication systems in terms of capacity, outage 
probability and symbol error probability (SEP) dramatically. 
Two main relaying protocols that have been researched a 
lot, are amplify-and-forward (AF) and decode-and-forward 
(DF) ID, El- In m-m performance of AF and DF protocols 
in terms of outage probability and SEP have been widely 
investigated over Nakagami-m and Rayleigh fading channels. 

In cooperative communication networks, the use of multiple 
relays to facilitate the source-destination communication was 
proposed to increase the spatial diversity gain Q, El- To avoid 
the interfering, source and all the relay transmissions must take 
place on orthogonal channels. Thus multiple relay cooperation 
is considered inefficient in terms of channel resources and 
bandwidth utilization. To overcome this problem, opportunistic 
relaying (OR) has been proposed at which only the best-relay 
from a set of L available relays is selected to participate in 
communication ED, GDI. It was shown that OR achieves full 
diversity order. With this technique, the selection strategy is to 
choose the relay with the best equivalent end-to-end channel 
gain which is obtained as the highest minimum of the channel 
gains of the first and the second hops under DF protocol or 
with the best harmonic mean of both channel gains under AF 
protocol fTR . lH2l . 

in m-m performance of Log-normal fading channels 
over different structures and relaying protocols has been 
investigated. From the practical point of view. Log-normal 
distribution is encountered in many communication scenarios. 
For instance, when indoor communication is used at which 
users are moving. Log-normal distribution not only models the 
moving objects, but also the reflection of the bodies. Moreover, 
it models the action of communicating with robots in a closed 


environment like a factory ED. In indoor radio propagation 
environments, terminals with low mobility have to rely on 
macroscopic diversity to overcome the shadowing from indoor 
obstacles and moving human bodies. Indeed, in such slowly 
varying channels, the small-scale and large-scale effects tend 
to get mixed. In this case. Log-normal statistics accurately 
describe the distribution of the channel path gain El. 

A main underlying assumption in majority of the current 
literature on cooperative communication is the availability of 
the channel state information (CSI) at the receiver. Recently 
there has been an interest in evaluation the performance of re¬ 
lay networks over imperfect channels El-El. Performance 
analysis of opportunistic AF and DF schemes over Log-normal 
channels is a non-trivial task when the CSI is imperfectly 
known at all nodes. To the best of our knowledge there have 
been no reported results on the outage probability analysis of 
such a systems yet. The main contribution of this paper is 
to derive closed-form expressions for the outage probability 
of the multi-relay DF and AF systems, employing “best- 
worst” and “best harmonic mean” relay selection criteria over 
imperfect Log-normal fading channels, respectively. 

The rest of the paper is organized as follows: Section HU 
describes the system model. SectionQjl] discusses the Gaussian 
error model for imperfect channel estimation. We take advan¬ 
tage of this model in our relay selection scenario to obtain 
output instantaneous SNR. Section [Iv] provides the harmonic 
mean of two Log-normal random variables (R.V) for AF and 
the equivalent CDF of the best-worse selection criterion for 
DF in order to calculate outage probability. SectionlVlpresents 
simulation results, while Section IVll provides some concluding 
remarks. 

II. System Model 

We consider a multi-relay scenario, in which a source node 
(S) communicates to a destination node (D) via multiple fixed 
relays ( Ri , l = 1,L). We assume that there is no direct link 
between the source and the destination, and communication 
occurs using a two-hop protocol over two time slots El- 
The fading coefficient over source to relay l, and relay l 
to destination are denoted with h s , ri and h rit d which are 
assumed to be independent and non-identically distributed 
Log-normal R.Vs. Dropping the indexes, channel gains during 
the transmission of a bit are modeled by ft, = 10 o lJf , at which 
X ~ A'T/J, (T 2 ) and Aft denotes a Gaussian distribution. During 
the first time slot, source broadcasts the signal to L relays. In 


the second time slot, only the best-relay forwards the signal 
to the destination, and source remains idle. Let us denote with 
P s the power transmitted by the source and thus the set of 
below equations summarize the operation taking place for each 
symbol 

JJs,ri — T tl ri l lj * * * j L (1) 

Vri,d — %rh'r[,d A tldi (2) 

where x is the transmitted signal with power .E[|x| 2 ] = 1, 
and n r j and n,i are complex additive white Gaussian noise 
(AWGN) in the relay and destination, respectively. Without 
loss of generality, we assume that all the AWGN terms have 
equal variance as Nq. For AF relaying x r = Aiy sri at 
which relay amplification factor, Ai, is chosen to satisfy an 
average power constraint and will be defined later. For DF 
relaying x r = \fT\x, where x is obtained after demodulating 
y s ,n followed by modulating for retransmit to the destination. 
Therefore, in brief we study the OR with two conventional 
relaying strategies at the relay: 

• DF: The best-relay decodes the message, re-encodes it 
and transmits that message in next time slot. 

• AF : In the second time slot the best-relay process the 
received signal and forwards it to the destination. 

In the sequel we investigate the performance of OR-DF 
and OR-AF schemes over Log-normal fading channels where 
relay and destination are provided with a estimation of their 
corresponding channels. 

A. Transmission with DF Protocol 

In the second time slot, for DF scheme, only the best-relay 
according to the best-worse criterion El is chosen to decode 
and re-encode the received signal which yields x. Then the 
selected I th relay send x r = \/T\ x to the destination where 
P r is the relay power. As a result the received signal at the 
destination is given by 

Vri,d — sjPrXh ri _rl T Tldi (3) 

B. Transmission with AF Protocol 

Under AF protocol the relay with the best harmonic mean 
of both source to relay and relay to destination gains is 
chosen to forward x r = Aiy S)ri to destination 0. Note that 
Ai = yj ^ asserts the relay amplification factor 

which controls the output power of the relay (TJ. Since A; 
depends on the fading coefficient, each relay has to estimate 
its own received channel. We assume that relays estimate their 
corresponding channels, h Stri , and then use it to amplify the 
received signal. The received signal at the destination is of the 
form 


III. Instantaneous SNR with Imperfect CSI 
We denote the estimated and exact channel coefficients as 
h and h, respectively. To estimate the h linearly with respect 
to h, we employ the following model l 5Tft 

h = ph + e, (5) 

where e is the channel estimation error modeled by zero mean 
complex Gaussian distribution with variance cr 2 and p is the 

correlation coefficient between h and h which is given by 

2 

p = £ ED- The variances of the error and the exact channel 

a h 

coefficient are related by p as 

= ( 6 ) 

Since estimated channel is the combination of a Log-normal 
and a complex Gaussian R.V, it’s pdf does not exactly follow a 
Log-normal distribution. However, in Section [V] we will show 
that when the correlation coefficient between the estimated 
channel and the real one is near to one ( i.e ., p ~ 1), 
this approximation is acceptable. By employing the imperfect 
channel estimations at the receiving node we obtain a closed- 
form expressions for instantaneous SNR at the destination. 


A. Transmission with DF Protocol 

Substituting (0 into 0 and 0, we have 

Vs,ri—sj l *S Xf> s r: h St ri+ \jP$XG S ^riA tin l — 1; L (7) 
Uri,d — -ttrPri,dhs 7 ri "F X r C ri ^d A Tld■ (8) 

After receiving the signal in the I th relay, since the noise power 
is not the same on all sub-channels, each diversity branch 
has to be weighted by its corresponding complex fading gain 
over total noise power on that particular branch. Therefore, 
the selected relay and the destination will decode the received 
signal using MRC as 0 


h* h* 

rl s,r, - n n,d 

Vs,ri — Ar Us,rn Vri,d — AT Uri,dy 

1\ 0 


(9) 


then using 0, instantaneous SNR at the relay and destination 
can be formulated as 

p „ 

Isf/i = T s,ri-jrr\h S: r i\ = Ts,n7s,ri |^s,r; | = T s, ri Ts,rn (10) 

jV o 


Trf.d~ T n,d AT \hri,d\ — T ri ,dHn ,d\hn ,d\ — Tri,dTri,d, (H) 

An 


Pr 


i ^s,rj \ IL s,ri | , Pri,d\hri,d\ 

where r s r , = ——- - — or,H " r - — — -— 


1 s,n — /■, 0 j ano x ri d — 


i+- 


N 0 


1 +- 


B. Transmission with AF Protocol 


Remembering that our transmission model for AF scheme 
is given by 0, substituting 0 into 0, the received signal in 
the destination is 


Uri ,d — Aiy Stri hr u d T tld 

— Al h s ^n fori , d-t' T Aih ri ’dtlri A tld- (4) 

Armed with these system models, in the consecutive sections 
we model the imperfect CSI at the receiving nodes, and then 
study the performance of aforementioned AF and DF schemes 
with relay selection, in term of outage probability. 


Vri,d 


y^Aix(ps, ri Pri ,dhs,ri hri ,d) 

's/Ps■ / ^-lX^Ps,n^ j s,ri^s,ri H - Pri^dJ^ri^d^ri^d H - &s,r 
■AlPn : dhri,dCri,d H - Pri,d^ri ,d^r H - Tld 5 


&ri ,d)~^~ 

( 12 ) 


where, the first term presents the received signal, the second 
and third terms stand for the error signal, and the overall noise 
















at the destination is h d — Aip ru dh ri ,dn r + ri d which is a 
complex Gaussian R.V with fid ~ A/”(0, cr|) where a? is 


a? = Nn 



Pr Pri ,d | h'ri ,d \ 
Ps\k,r, | 2 + N 0 


(13) 


Using MRC at the input of destination, the estimated signal is 


Vri,d 


h* h* 

L s,ri l ri,d 


Dn,d- 


(14) 


Supposing that n r , rid, e Sjr , and e ri , d are processes that 
are independent from each other, the instantaneous SNR at 
the destination is obtained as o. at the top of the next 
page. In order to have a more tractable form, we neglect 
N o(^WE a e a r a e r d + 1) i * 1 <n~5l . So, we can further simplify 
02 ) as 


- eff 

id,i 


^s,ri'fs,ri 4“ ^ri,d^)ri,d 


(16) 


where A Sin = 1 


P,\h 


I 2 - Pr\h. 

N 0 > 7 n ,d~ 


N 0 


l,d, X rit d — 1+ P rt ,d e s,ri, 7s,r; — 

r Pr^ 2 , 

’ £ «.n = N ’ and e r u d=- Nq 


e ri,d 


IV. Performance Analysis 


In this section we derive closed-form expressions for outage 
probability of the OR scheme under AF and DF protocols. 
Outage probability is defined as the probability that the instan¬ 
taneous SNR at the receiver, 7, falls below a predetermined 
protection ratio, 7 th , namely 

[1th 

Pout = P [7 < 1th] = / fr(l)di. (17) 

Jo 

where represent the pdf of the instantaneous SNR. It 

can readily be seen that the outage probability is actually the 
cumulative distribution function (CDF) of 7 evaluated at 7 th- 
Before proceeding, we introduce following theorems from l22l 
which will be used in the sequel to derive /x( 7 ). 

Theorem 1: If X and Y are two R.Vs with relation Y = 
m.X, then /y( 7 ) = F fx(%). 

Theorem 2: If X is a Log-normal R.V with distribution X ~ 
LogN(p,x,cr‘x), then Y = i is a Log-normal R.V distributed 

as LogN(-px,Vx)- 

Theorem 3: If X is a Log-normal R.V with distribution 
X ~ LogNijix ■ p’x), then Y = X 2 is a Log-normal R.V with 
distribution defined as f~ LogN(2px,^<r\). 

Theorem 4: If X is a Log-normal R.V with distribution X ~ 
LogN(px,crj c ), then considering theorem 1, Y = mX is a 
Log-normal R.V with distribution defined as Y ~ LagN (jjx + 
10 log {m),<r%). 


A. Transmission with DF Protocol 

If we suppose that h s>ri and h ri}d , the estimations of the 
source-relay and relay-destination channels, respectively are 
available at destination, in the first phase, destination node 
equipped with selection combiner (SC), selects the worst hop 
of each branch as 

leqi = min {7.5. r( , X fr l: d}- (18) 


In second phase, SC selects the branch with the best 7 eqi as 

7 SC = max{7 egi , % q2 ,..., % qL }. ( 19 ) 


Since the SC chooses the weakest part of each branch and 
then the best one is selected to send the signal, the occurrence 
of outage is equal to the case when the best weak link’s SNR 
is under the threshold (7 t h). 

Pout = P(.Tsc < 7 th)- ( 20 ) 


Without loss of generality, we stipulate equal power alloca¬ 
tion to source and best-relay (P s = P r = P). By considering 
the pdf of Log-normal R.V (23], h, with corresponding Normal 
parameters defined as //.^ and a?, then by using theorems 3 
and 4, respectively, after some elementary manipulations and 


dropping the indexes, the pdf of 7 = 
obtained as 


_p 

N 0 


= 7 


ft (7) = 


£ 


■ exp 


(101og 10 7-Pr) 2 


Mf = 10 iogio-^TO-S log 10 tf (7), 


2 a? 

T 

100 

^=hns 4,(T) - 


( 21 ) 

( 22 ) 


where, £ = = 2p~ h + 10 log 7 , a? = 4a?, 

^( 7 ) = (l + V ariff) / E(p/) 2 }. Relations (l22l > help us derive 
the parameters of Log-normal distribution directly from the 
variable, 7. We also can express the CDF of 7 as C3 


J 1 f(7) = g( Mt 1 g 01 ° SlO7 )> 7>0 (23) 

2 

where, Q(x) = -^= e~ du is the standard one¬ 

dimensional Gaussian function. According to independency of 
channel gains, the pdf of 7 eqi in (fl 8 l > can be expressed as 


Pf eqi (7) = P (min{7s,n, 7 n,d < 7» 

= 1 -( 1 - F r a , n (W( l - F t ri , d m- (24) 

Substituting (1231) into ( 1241 ). yields 


F f (7) = 1 — (1 — Q (f2 s ri ))(l — Q (Q ri d )). (25) 

e< ll 

where D ai b = ([if ^ — 101 og 10 7 )/a 7 ( for a £ {s,n} and 
b £ {ri,d} . According to the fact that our branches are 
independent and using ( fl9l >. P out is given by 


Pout — P(,Teqi ’Fi^ith)P{f)eq ‘2 — Tth') * * * P(,1feqL — Tth)- (26) 

Substituting (1251) into (l26l > yields outage probability for OR- 
DF scheme as 

L 

Pout — n (Q(n s>ri )+Q(n ri>d )-Q(n s , ri )Q(n ritd )). (27) 

1=1 

B. Transmission with AF Protocol 

In the OR-AF scenario, destination selects the maximum 
harmonic mean of both S-R and R-D channel gains. As a 
result the outage probability can be readily obtained as 

Pout = P (lsc= max{ 7 ^/, 7 e JJ , • • • , 7 ^/} < 7 th) , (28) 

Similar to DF mode, using the independency of the branches, 
we arrive at 

Pout = P{l e d,I<nth)p{l e d f J <7 th) • • -P(p/doL^yth) , (29) 



















(15) 


~ e ff 

ry — 

id,l ~ 


Ps PrPs,riPri,d\hri,d\ |^s,r' 


^0 2 (fl^,nl 2 (l + p: 


s,ri No 


)<,« + 


Pr. 

N 0 


I K,d\ 2 (l 


o 2 ^ 
Pr L ,d N 0 


K 2 


N 0 


+ 1 

r l Mq e ri ,d 


which is the product of CDFs of ■ l = 1, • • ■ , L. The CDF 
of is given in following preposition. 

Proposition 1: The pdf of the received instantaneous SNR 
at the destination for OR-AF relaying protocol over imper¬ 
fect non-identical Log-normal fading channels is %/n ~ 

LogN {-p x + 10log P CSr^l a ^- 

Proof: We can rewrite (Il6l > as 

-eff P s,r\P ri ,d 1 

7 d ,A = 1 - \ - \ - Z — + T- Z - 

In order to simplify the analysis we introduce a new random 
variable \ given by the summation of two random variables as 

y = a + B, where a = -r—4- and B = -r—- . Remem- 

bering that the pdf of the S-R and R-D channels are given by 
hs, n ~ LogN (p s ,ri, cr 2 r ,) and h n , d ~ LogN (p n ,d,c r* jd ), 
it become obvious that both a and f} have a Log-normal dis¬ 
tribution as a ~ LogN (-(2// Siri + lOlog(%, n A s , r J),40f ri ) 
and (3 ~ LogN (~(2p ri:d + 10 log( 7 rijd A njd )), 4cr^.d)- Next 
we derive the pdf of For this purpose we employ the 
Wilkinson method which is described in m as follows: 

Wilkinson method : If X\, ■ • • , Xn are N Log-normal R.Vs, 
then Z = X; + X 2 + • • • + A7\-, can be approximated by a 
Log-normal R.V with the following parameters 

= i ^21n(«i)-i ln(u 2 )Y oz = ^ v/ln(M 2 ) -2 ln(wi) (31) 



u± = ^2 e^ Xi + ' 


L ) _ e Wz+^f] 


(32) 


i=l 

N 

u 2 = Y J e 

i =1 


N—l N 

(2px i +2crx i ) _|_ 2 ^ ^ ^ ^ 

2=1 _ 7 = 2 +l 




+2 r ij <T x i <T x j ) _ e (2fL Z +2(J 2 z ) 


(33) 


where A = ^ 7 ^ and rij is the correlation coefficient between 
Xj and X , which is defined as 


E{(Xi - pxi){Xj - p Xj )} 


VXi&Xi 


(34) 


Replacing p Xl , <j Xl , p X2 and cr X2 with p a , cr a , pp and 
<jp and using ( Tut . p x and o' 2 are obtained. Finally, using 
theorem 2 and 1, respectively. It can be readily checked that 
7 has a Log-normal distribution as ~ LogN {—p x + 


10 log ■ 


-■(J 2 ). As N becomes large, the central limit 


theorem states that the sum will become close to Gaussian 
distribution. ■ 

Using (l23l ). the outage probability of the OR-AF relaying 
over imperfect non-identical fading channels can be expressed 
as 


L p^tf - 10 log 10 7 

Pout = Y[Q(^L_ -). 


1=1 


(J - eff 

Id’n 


(35) 



Fig. 1. Outage probability for OR-DF mode over imperfect Log-normal 
fading channels for different p. 



Fig. 2. Outage probability for OR-AF mode over imperfect Log-normal 
fading channels for different p. 

V. Numerical Results 

In order to justify our analytical results, we provide some 
monte-carlo simulations in this section. In all simulations, 
unless it mentioned otherwise, we have three relays {L = 3), 
7 th = 3, and p s>ri = p rt ,d ■ Log-normal fading channel 
parameters and the correlation coefficient between the exact 
and estimated channel, p, are given in each figure. We can 
calculate the variance of the each channel (o 2 ) using l22ll 

a 2 =p 2 x ( lO^too-l), (36) 

where, p x and n 2 x are the Log-normal fading parameters. 

Fig HI and Fig. [ 2 ] depict the effect of channel estimation 
error on OR-DF and OR-AF protocol over Log-normal fading 
channel, respectively. 07 1; pip, <ji 2 , pi 2 are the parameters of 
the first and second hop, respectively. Following conclusions 
are drawn from Figs. HI and [ 2 ]: 

1) As the p decreases, or equivalently, the estimation error 
increases, the performance become worse. 

2) As p decreases the distance between the simulation and 
analytical result increases. 

3) Depending on p, in a specific SNR, the performance of 
the system saturates, that is by increasing SNR, we do 
not get improvement in the system performance. 

The second conclusion, is the outcome of the approximation 
that was mentioned in section Ell However, simulations show 
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Fig. 3. Comparing the accuracy of the analytical outage probability derived 
for OR-AF using Wilkinson method with the monte-carlo simulation. 



Fig. 4. Comparison between analytical and simulated outage probability 
of OR-DF and OR-AF scheme with imperfect CSI with different number of 
relays. (L = 2, 3) 


that approximating the estimated channel as a Log-normal R.V 
is acceptable. We observe that for p ~ 1, the analytical result 
is acceptable for good range of SNRs. 

Fig. [3] investigate the accuracy of Wilkinson method with 
different variances for OR-AF system. As it is shown in ll24l . 
we see that the higher the variance, the higher the difference 
between the parameters of the approximated pdf according to 
Wilkinson method and the pdf according to simulation. 

As it has been mentioned in J4), we can see that in Fig. [4]the 
DF outperforms the AF protocol in low SNRs; however, they 
become close to each other in high SNRs. We can also notice 
that the number of cooperating relays has a strong impact on 
the performance enhancement, that is the achieved diversity 
order is related to relay number, L. 

VI. Conclusion 

In this paper, end-to-end outage probability of a wireless 
communication system using best-worse and best-harmonic 
mean selection Over Imperfect Non-identical Log-normal Fad¬ 
ing Channels has been investigated. Channel estimation error 
has been considered as a Gaussian random variable. Since the 
distribution of estimated channels are not Log-normal either 
as would be in the case of the Rayleigh fading channels, 
simulation results do not exactly follow the analytical results. 
However, this difference is negligible for a wide variety 
of situations. Harmonic mean of two Log-normal R.Vs in 
presence of channel estimation errors was also derived. 
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